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ABSTRACT 

Based on the recent very deep near-infrared imaging of the Hubble Ultra Deep Field with WFC3 on the 
Hubble Space Telescope, five groups published most probable samples of galaxies at z ~ 8, selected by the so- 
called dropout method or photometric redshift; e.g., F^-dropouts (fios - ./125 > 0-8)- These studies are highly 
useful for investigating both the early star formation history of galaxies and the sources of cosmic re-ionization. 
In order to better understand these issues, we carefully examine if there are low-z interlopers in the samples of 
Z ~ 8 galaxy candidates. We focus on the strong emission-line galaxies at z ~ 2 in this paper. Such galaxies 
may be selected as Fio5-dropouts since the [Olll] A 5007 emission line is redshifted into the /125-band. We 
have found that the contamination from such low-z interlopers is negligibly small. Therefore, all objects found 
by the five groups are free from this type of contamination. However, it remains difficult to extract real z ~ 8 
galaxies because all the sources are very faint and the different groups have found different candidates. With 
this in mind, we construct a robust sample of eight galaxies at z ~ 8 from the objects found by the five groups: 
each of these eight objects has been selected by at least two groups. Using this sample, we discuss their UV 
continuum slope. We also discuss the escape fraction of ionizing photons adopting various metallicities. Our 
analysis suggests that massive stars forming in low-metallicity gas (Z~5x IO^Zq) can be responsible for the 
completion of cosmic re-ionization if the escape fraction of ionizing continuum from galaxies is as large as 0.5, 
and this is consistent with the observed blue UV continua. 
Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

In order to understand both the early phase of galaxy evo- 
lution and the cosmic history of intergalactic space (i.e., the 
cosmic re-ionization of intergalactic medium), it is crucially 
important to probe any objects at very high redshift. In this 
decade, ground-based 8-m class telescopes have enabled us 
to find a large number of galaxies beyond z ~ 6 (e.g., Hu et al. 
2002; Rhoads et al. 2004; Taniguchi et al. 2005; Stern et al. 
2005; Kashikawa et al. 2006; lye et al. 2006; Ota et al. 2008; 
Ouchi et al. 2009; see for a review Taniguchi 2008). Also, 
the optical deep imaging survey promoted with the Hubble 
Space Telescope (HST), the Hubble Ultra Deep Field (HUDF; 
Beckwith et al. 2006), resulted in very probable detection of 
galaxies beyond z = 7 together with deep near-infrared (NIR) 
imaging surveys with HST and large ground-based telescope 
facilities (Bouwens et al. 2006, 2008). In addition, recent 
deep near infrared imaging surveys with 8-m class ground- 
based telescopes provided probable candidates of galaxies be- 
yond z = 7 (e.g., Hibon et al. 2010). 

More recently, the successful installation of the Wide Field 
Camera 3 (WFC3) on HST was promptly used to observe 
HUDF in the NIR window. Using these imaging data, very 
promising candidates of high-z galaxies at z ~ 7 - 10 were 
found by several researchers (Bouwens et al. 2009, 2010b; 
Oesch et al. 2010; Bunker et al. 2010; McLure et al. 2010; 
Yan et al. 2010; Finkelstein et al. 2010). These objects 
were used to investigate both the early star formation his- 
tory in galaxies and the cosmic re-ionization history. In fact, 
the following new results were obtained in these studies. (1) 
The number density of luminous galaxies appears to become 
smaller with increasing redshift while that of less luminous 
ones appears to be constant (Bouwens et al. 2010b). (2) 
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The galaxies beyond z ~ 7 tend to have bluer UV continuum 
with (3 < -3 (f\ oc A' 3 ), significantly bluer than the typical 
/3 ~ — 2 of nearby starburst galaxies and Lyman break galax- 
ies (LBGs) at z < 6 (Bouwens et al. 2010a, 2010b; Finkelstein 
et al. 2010). This trend is more pronounced for fainter galax- 
ies. And, (3) the escape fraction of ionizing UV continuum is 
also discussed in terms of the cosmic re-ionization at z > 7. 
For example, if more numerous, low-luminosity star-forming 
galaxies could be present, they would be enough to reionize 
the intergalactic medium together with both a larger escape 
fraction of ionizing photons from them and a small clump- 
ing factor of the universe (e.g., Yan et al. 2010; Bunker et 
al. 2010). However, there appears no consensus among the 
studies. 

Any conclusions and suggestions derived through these 
studies depend strongly on the robustness of the detection of 
such very high-z galaxies. In order to examine how the results 
are robust, it is important to reject any contamination from 
low-z interlopers. Here, we focus on galaxies at z ~ 8 (i.e., 
F105 dropouts) because the following five groups extracted 
their own samples of galaxies at z ~ 8 based on the WFC3 
observations of the HUDF; Bouwens et al. (2010b), Bunker 
et al. (2010), McLure et al. (2010), Yan et al. (2010), and 
Finkelstein et al. (2010). 

Possible interlopers are low-temperature stars, low-z red 
galaxies, and low-z very strong emission-line galaxies with 
faint continua. All groups but Bunker et al. (2010) removed 
compact objects (i.e., possible low-temperature stars in our 
Galaxy) from their samples. They also checked the possibil- 
ity of low-z red galaxies by using Spitzer IRAC photometry 
and confirmed that such contamination cannot occur in their 
analysis. However, they did not examine the possibility of 
contamination from strong emission-line galaxies at low z, al- 
though Oesch et al. (2010) pointed out such possibility only 
for z ~ 7 (zg5o -dropout) galaxies. Therefore, in this paper, we 
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first examine whether or not this type of contamination occurs 
in their analysis. 

Then, we construct a robust sample of galaxies at z ~ 8 
free from any contamination. In this procedure, we also take 
account of possible false detection due to unexpected photo- 
metric effects discussed by Bouwens et al. (2009). Based on 
this new sample together with the results of Bouwens et al. 
(2010b, 2010c), we discuss the early star formation history of 
galaxies and cosmic re-ionization at z ~ 8. 

Throughout this paper, magnitudes are given in the AB sys- 
tem. We adopt a flat universe with £l ma tter = 0.3, = 0.7, and 
//o = 70kms" 1 Mpc" 1 . 

2. THE SELECTION CRITERIA FOR IDENTIFYING GALAXIES AT 

Z~8 

Before studying the possibility of contamination by strong 
emission-line galaxies at z ~ 1 .5 - 2, it is necessary to exam- 
ine how the galaxies at z ~ 8 have been identified in HUDF 
by the following five groups: (1) Bouwens et al. (2010b), (2) 
Bunker et al. (2010), (3) McLure et al. (2010), (4) Yan et al. 
(2010), and (5) Finkelstein et al. (2010). In this section, we 
summarize their selection criteria for identifying galaxies at 
Z ~ 8. 

First, we list up all the candidates identified by the five 
groups in Table 1 . The numbers of detected galaxies at z ~ 8 
are 5, 7, 8, 15, and 9 in Bouwens et al. (2010b), Bunker 
et al. (2010), McLure et al. (2010), Yan et al. (2010), and 
Finkelstein et al. (2010), respectively. Their studies lead to 
20 galaxies at z ^ 8 in total. Note that the last object (ID 
No. 21) is selected as a zsso dropout in Oesch et al. (2010), 
Bunker et al. (2010), and Yan et al. (2010) while as a z ~ 8 
galaxy candidate in McLure et al. (2010) and Finkelstein et 
al. (2010). We treat this object as a galaxy at z ~ 7 and thus 
we do not use it in this paper. Yan et al. (2010) identified their 
z8-SD15 as YD5 in Bunker et al. (2010). However, checking 
their positions and thumbnails, we identify z8-SD24 as YD5 
and there is no counterpart of z8-SD15 in the samples of the 
other three groups. 

Only four galaxies (Nos. 1, 2, 3, and 4 in Table 1) are 
commonly identified by all five groups. Since the five groups 
made their own data reduction, their final images are different 
from each other. In fact, their zero points are slightly differ- 
ent from each other. In addition, their selection criteria for 
galaxies at z ~ 8 are also different; see a summary of their 
selection criteria given in Table [2] Three groups (Bouwens et 
al. 2010b; Bunker et al. 2010; Yan et al. 2010) selected z ~ 8 
galaxies using the dropout method (Fi 05 -dropout) while the 
others (McLure et al. 2010; Finkelstein et al. 2010) selected 
z ~ 8 galaxies using the photometric redshift technique. Al- 
though the selection criteria are slightly different from group 
to group, we adopt F105 — /125 > 0.8 as our primary criterion 
for galaxies at z ~ 8. 

3. STRONG EMISSION-LINE GALAXIES AS INTERLOPERS 

Possible interlopers are very strong emission-line galaxies 
with very weak continuum emission like blue compact dwarf 
galaxies; e.g., I Zw 18 (Izotov et al. 1999). Strong emission 
lines may modify observed broadband colors and thus make 
false continuum break featurefl 

3 Here, we do not consider active galactic nuclei (AGNs). Although AGNs 
are also strong emission-line sources, they have brighter rest-frame blue con- 
tinuum and thus they may be easily detected in the very deep ACS imaging 
of HUDF. Therefore, we focus on star-forming galaxies like blue compact 
dwarf galaxies in the local universe. 



Actually, in the local universe, metal-poor galaxies with 
large emission-line equivalent widths have been surveyed us- 
ing their unusual broadband colors caused by strong emission 
lines (Brown et al. 2008). Recently, using the data of the 
Sloan Digital Sky Survey Data Release 7 (Abazajian et al. 
2009), Cardamone et al. (2009) investigated a class of com- 
pact star-forming galaxies called "Green Peas." Because of 
their very strong [Olll]A5007 emission line with a rest-frame 
equivalent width of EW ([Olll])> 500 A, "Green Peas" in the 
redshift range of 0. 1 12 < z < 0.360 are very bright in r band. 
Their r— i colors are bluer by > 0.5 mag than those of normal 
galaxies. Such ultra strong emission-line galaxies have also 
been found up to z ~ 1 (Kakazu et al. 2007; Hu et al. 2009; 
see also Ohyama et al. 1999). 

Since the strongest emission line is generally [O lll]A5007 
for such dwarf starburst galaxies (e.g., Izotov et al. 1999; 
Brown et al. 2008; Cardamone et al. 2009), it is possible that 
strong emission-line galaxies at z ~ 1.5-2 could be detected 
as F105 dropouts if their optical continuum emission is very 
weak. To demonstrate this, we show a typical spectrum of a 
strong emission-line galaxy at z= 1.5 (Figure [TJ. Therefore, 
in this paper, we first examine if this type of contamination 
could accidentally occur in the identification of galaxies at 
z ^ 8 in HUDF. 

4. MODELS OF LOW-Z STAR-FORMING GALAXIES 

We make our model spectral energy distributions (SEDs) 
for star-forming galaxies. First, we generate model galaxy 
SEDs with nebular continuum emission using the population 
synthesis model, STARBURST99 (Leitherer et al. 1999), 
adopting the mode of constant star formation (SFR) rate to- 
gether with both Salpeter initial mass function (IMF; the 
lower and upper masses are 1 M Q and 100 M©, respectively) 
and the metallicity of Z = O.2Z = 0.004. These models corre- 
spond to nearby blue compact dwarf galaxies with very high 
SFRs and sub-solar metallicities (e.g., Z ~ 0.1Z Q ), such as I 
Zw 18 (Izotov etal. 1999). 

The nebular continuum is important in this analysis since 
it affects the global shape of the SEDs, especially for young 
starburst galaxies (Schaerer2003; Bouwens et al. 2010a). We 
calculate the luminosity of the nebular continuum (free-free, 
free-bound, and two photon emission for H and He) assuming 
that the escape fraction of the ionizing continuum is / esc = 
0. This value is consistent with observations of star-forming 
galaxies at z ~ 1 - 3 (e.g., Malkan et al. 2003; Siana et al. 
2007; Steidel et al. 2001); where we expect strong emission- 
line galaxies as possible interlopers. 

However, since nebular emission lines are not included in 
the above SEDs generated by STARBURST99, we add some 
emission lines to the SED generated above. Our procedures 
are as follows. First, we calculate the luminosity of H/3 emis- 
sion line, L(H/3), using the following relation (Leitherer & 
Heckman 1995): 

L(H/3) (erg s" 1 ) = 4.76 x lO" 13 ^ (s" 1 ), (1) 

where A^, yc is the ionizing photon production rate estimated in 
each SED model. Second, we estimate other emission lines 
(e.g., [O II] A 3727, [O III] A 5007, and so on), adopting 
the observed emission-line ratios of I Zw 18 (Izotov et al. 
1999), which is one of the most metal -poor galaxies known 
in the local universe. Then, we add all emission lines with 
F{\)/F(ftf3) > 0.004 from [O II] A 3727 to [Ar III] A 7135. 
Note that the emission-line equivalent widths tend to be larger 
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with decreasing metallicity because here we use the fixed 
emission-line ratios of I Zw 18. 

Given the star formation history, the equivalent widths of 
emission lines vary with the age of a galaxy because the 
emission-line luminosities are proportional to the production 
rate of ionizing photons. Figure |2] shows the evolution of 
equivalent width of H/3, [O III] A5007, and [O ll]A3727 in our 
model SEDs. In the next section, we examine the variation of 
observed colors as a function of redshift for the model galax- 
ies with ages of 1, 10, and 100 Myr. The values of EW (H/3) 
are 520, 160, and 70 A for galaxies with ages of 1, 10, and 
100 Myr, respectively. 

5. RESULTS 

First, we examine if some low-z strong emission-line galax- 
ies accidentally satisfy the selection criterion for galaxies at 
z ~ 8 galaxies; i.e., F105-./125 > 0.8. In Figure [3] we show 
the variation of F105 —/12s colors of the model galaxies with 
age of 1, 10, and 100 Myr described in Section 4 as a function 
of redshift (0 < z < 4). It is shown that the reddest F105 — /125 
color is achieved at 1.47 < z < 1.81 for the model with age 
of 1 Myr. In this case, the /125 magnitude is dominated by 
the four strong emission lines of [O lll]A5007, [O lll]A4959, 
H/3, and H7 (see Figure [TJ. However, the color is as red as 
F105— /125 — 0.4 because [O ll]A3727 contributes to the flux 
in F105. Another red peak with F105 — /125 — 0.35 is seen at 
z ~ 0.83. This feature is due to redshifted Ha emission in 
/125. A broad red peak around z ~ 2.5 is due to [O ll]A3727 
in 7n5, but the color is only F105 —J\2S ~ 0.1. Here we note 
that the [O II] emission is generally weaker than [O III] in 
metal-poor, star-forming galaxies (e.g., Nagao et al. 2006). 
For the models with ages of 10 Myr and 100 Myr, the reddest 
color appears to be Y\ 05-/125 ~ 0.2 - 0.25 at most around 
1.8. 

In summary, strong emission-line galaxies at 1 .47 < z < 
1.81 give rise to F105 —J\25 ~ 0.3 with the reddest color of 0.35 
at z ~ 1 .54. These colors are much bluer than the primary cri- 
terion of Fio5-7i25 > 0.8 (e.g., Bouwens et al. 2010b). How- 
ever, at this stage, we cannot conclude that there is no con- 
tamination by strong emission-line galaxies around z ~ 1.5 
because the following two points must be also taken into ac- 
count; (1) reddening, and (2) photometric errors in the WFC3 
observations of HUDF. 

As for the reddening, it is known that nearby blue com- 
pact dwarf galaxies are reddened in the range of < Ay < 1 
(Kong et al. 2003). However, we do not know how star- 
forming galaxies at z ~ 2 galaxies are typically reddened. 
Since they may be more massive than nearby blue compact 
dwarf galaxies, they may be more strongly reddened. For 
example, the extinction in star-forming galaxies at z ~ 2 se- 
lected as BX galaxies ranges from Ay ~ to 1.5 (Shapley et 
al. 2005). Taking account of this, for the safety, we investi- 
gate the three cases; Ay = 0, 1, and 2 mag together with the 
reddening curve for starburst galaxies obtained by Calzetti et 
al. (2000). To examine the reddening effect, first we use the 
F105— /125 versus /125 -//i60 diagram because this diagram is 
used to identify galaxies at z ~ 8 (e.g., Bouwens et al. 2010b). 
In addition, we also use another diagram of F105-/125 ver- 
sus 2775-/125 diagram. The reason is as follows. If some 
z ~ 8 galaxy candidates are strong emission-line galaxies at 
z ~ 2, the optical Advanced Camera for Surveys (ACS) imag- 
ing data of the HUDF are also useful. Since the ACS survey 
depth (5a, 0.35" 4>) is 29.4 AB for B435, 29.8 AB for V 606 , 
29.7 AB for i 775 , and 29.0AB for z 850 (Bouwens et al. 2010b), 



'775 — J125 gives the strongest constraint on the contamination 
possibility; see Appendix A. Following the selection criteria 
of z ~ 8 galaxies in the five groups, we calculate a lower limit 
of i'775 — /125 color adopting 2a upper limit of 2775. 

In Figure HI we show comparisons between our SED mod- 
els with the effect of reddening and the WFC3 photometry by 
the five groups. Since we also show the photometric errors 
in each diagram, we can examine the possibility of contami- 
nation from strong emission-line galaxies at z = 1.47 to 1.81 
taking account of both the effects of reddening and photomet- 
ric errors simultaneously. 

It is extremely difficult to estimate the properties of the pop- 
ulation of potential contaminant 1.5 < z < 2 emission line 
galaxies, but we can qualitatively estimate the robustness of 
each z ~ 8 candidate using the reddening required for an 
emission-line galaxy at 1.5 < z < 2 to mimic that object's 
colors. At first, we summarize our results as follows. (1) 
The z ~ 8 galaxy candidates identified both by Bouwens et 
al. (2010b) and by Bunker et al. (2010) are free from the 
contamination of strong emission-line galaxies unless strong 
emission-line galaxies with Ay > 2 exist at z ~ 1.5-2. (2) 
However, if strong emission-line galaxies with Ay = 1 - 2 are 
present at z ~ 1.5 - 2, it is possible that some of the z ~ 8 
galaxy candidates identified by McLure et al. (2010), Yan et 
al. (2010), and Finkelstein et al. (2010) are strong emission- 
line galaxies at z ~ 1.5-2. 

Here we should note that the above different results among 
the five groups are due to the differences in photometry; i.e., 
some candidates have a non-zero probability to be an inter- 
loper in the sample of one group, but not another. How- 
ever, more importantly, we should also mention that heavily 
obscured dwarf galaxies with very strong emission lines are 
probably rare from observational material obtained to date. 
Therefore, we conclude that the possibilities of low-z inter- 
lopers are negligibly small for all the z ~ 8 galaxies identified 
by each group. 

6. DISCUSSION 

6.1. A sample of robust galaxies at z ~ 8 

In this section, we discuss the properties of galaxies ro- 
bustly determined to be at z ~ 8. In the last section, we 
have shown that the contamination of low-z strong emission- 
line galaxies is negligibly small for all candidates studied by 
the five groups. However, it is important to note that there 
might be other unexpected contamination sources (see, e.g., 
the summary of Bouwens et al. 2009), and so careful proce- 
dures are absolutely necessary to identify robust candidates at 
high redshift even if we use deep WFC3 imaging data. In par- 
ticular, they find that a large fraction of nine candidates (10 - 
18 in Table [U identified only by Yan et al. (2010) are quite 
close to foreground galaxies and they may be affected by fore- 
ground bright sources. We therefore do not use the nine ob- 
jects. We also do not use both 19 and 20 that are detected only 
by Finkelstein et al. (2010) since these objects are also close 
to foreground bright objects. Taking all these into account, 
we have constructed a sample of eight robust galaxies at z ~ 8 
all that are identified by at least two groups. These galaxies 
are from 1 to 8 in Table 1 . We summarize their photometric 
properties in Table 3. We note here that our robust sample 
is merely robust in the sense that these candidates have been 
selected by at least two groups and do not lie near foreground 
galaxies: spectroscopic follow-up is necessary to confirm that 
they are truly robust z > 8 galaxies. 
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Here it is noted that all the five groups tried to do their 
best in extraction of the total magnitude for their sample 
galaxies although their methods and photometric zero points 
are slightly differenfl In fact, the magnitudes and colors 
for the same objects tabulated in Table 3 are slightly dif- 
ferent among the groups. Therefore, in our later analysis 
(Section 6.3), we use mean values given in the last column 
of Table 3. For the value of Xzt a, we calculate the aver- 
age of x ((x)), the systematic error of x (<7 sys ), the mean of 
the random error a r d m , and the total errors <r tot are as fol- 

lows: (x) = (Em xd/N, a sys = ^Jd/N)^^- (x))\ a tim = 

\/ (Ef=i °f )/#> and o'tot = y^sys + ofdm' where N is the num- 
ber of data. The lower limit values are simply ignored for 
calculating the mean photometry, although in all cases they 
are consistent with the mean value from the true detection. 

6.2. Comment on the cosmic re-ionization 

In this section, we discuss how the star-forming galaxies 
contribute to the cosmic re-ionization at z ~ 8. Since the cos- 
mic reioniztion is one of important issues in high-redshift uni- 
verse, several investigations have been made for these years. 
Based on the survey of z-band dropouts in the Subaru Deep 
Field and the Subaru XMM-Newton Deep Survey, Ouchi et 
al. (2009) evaluated the ionizing photon production rate of 
star-forming galaxies at z ~ 7 adopting Salpeter stellar initial 
mass function (IMF) and sub-solar metallicity (Z = 0.2Z Q ). 
Comparing it with the critical ionizing photon production rate 
obtained by Madau et al. (1999), they found the ionizing 
photon production rate of star-forming galaxies is enough to 
complete re-ionization of the universe if the escape fraction is 
larger than 0.2 and suggested that the properties of very high 
redshift (z ~ 7) galaxies are different from those at low red- 
shifts, e.g., a larger escape fraction, a lower metallicity, and/or 
a flatter initial mass function if the very high redshift universe 
is ionized by only galaxies. Bunker et al. (2010) and McLure 
et al. (2010) also discussed the relation between the high-z 
galaxies and the cosmic re-ionization at z ~ 7 based on the 
survey of z-band dropouts in the Hubble Ultra Deep Field. 
Bunker et al. (2010) pointed out that the ionizing photon pro- 
duction rate of star-forming galaxies at z ~ 7 is lower than the 
critical value adopting the relation for the solar metallicity. 
McLure et al. (2010) calculated the evolution of the filling 
factor of the ionized hydrogen assuming the sub-solar metal- 
licity (Z = 0.2Z©) and the escape fraction of 0.2 and found 
that the inter-galactic medium would not achieve a unity HII 
filling factor until z = 4.2. Their findings also suggested that 
the star-forming galaxies at very high-z are the large escape 
fraction of ionizing photons, the lower metallicity, and/or the 
top-heavy IMF. 

Now, we focus on the cosmic re-ionization by galaxies at 
z ~ 8. Here we use the formalism studied by Madau et al. 
(1999). However, we do not use the critical SFR to complete 
the cosmic re-ionization (e.g., Equation (27) in Madau et al. 

4 The photometry methods of the five groups are summarized below; (1) 
Bouwens et al. (2010b): MAG_AUTO with Kron factor = 1.2 and then cor- 
rected to the case of Kron factor of 2.5 (this corresponds to the so-called 
aperture correction), (2) Bunker et al. (2010): aperture photometry with 0.6 
arcsec <f> and then applied an aperture correction, (3) McLure et al. (2010): 
aperture photometry with 0.6 arcsec (j> and then applied an aperture correc- 
tion, (4) Yan et al. (2010): MAG_AUTO with Kron factor = 1.2 and then ap- 
plied an aperture correction, and (5) Finkelstein et al. (2010): MAG_AUTO 
with Kron factor =1.2 and then applied an aperture correction. 



1999) since the Salpeter IMF and the solar metallicity are as- 
sumed there. Instead, we estimate the production rate of ion- 
izing photons based on the UV luminosity functions obtained 
by Bouwens et al. (2010b, 2010c). Using Equation (26) in 
Madau et al. (1999) together with the update of cosmological 
parameters given in Section 1, we obtain the following critical 
ionizing photon production rate density, 

• _ 1.2 x lO^photonss" 1 Mpc" 3 /l+z\ 3 / Vl b h 2 1Q \ 2 fC\ 

10n " U \~) 1004637 [sj 

(2) 

where is the cosmic baryon density and /770 is the Hub- 
ble parameter in units of h = 0.7. C is the clumping factor 
of neutral hydrogen, C= (Phi)(phi)~ 2 , where pm and (pm) 
are the local and cosmic mean density of neutral hydrogen 
(Gnedin & Ostriker 1997). We adopt Q h h% = 0.0463 from 
Wilkinson Microwave Anisotropy Probe 5 Year data (Komatsu 
et al. 2009). The clumping factor C is conventionally assumed 
to be 30 based on cosmological radiative transfer simulations 
by Gnedin & Ostriker (1997). However, smaller values for C 
have been recently suggested by cosmological hydrodynamic 
simulation in which photo-ionization heating together with 
radiative cooling is taken into account (Pawlik et al. 2009). 
Since the pressure support increases and then smooths out 
small-scale density fluctuations, the photo-ionization heating 
reduces the clumping factor down to C ~ 3 - 5. The clumping 
factor depends on the starting redshift of reheating, z r - If z r = 
19.5, C ~ 3, while if z r = 9, C ~ 5. In this paper, we adopt 
C = 5. Therefore, the production rate of ionizing photons is 
determined by only one parameter, / esc , that is the average es- 
cape fraction of hydrogen-ionizing photons from galaxies to 
intergalactic space. 

In Figure. [5] we compare the critical ionizing photon pro- 
duction rate density with the cumulative ionizing photon pro- 
duction rate densities corresponding to the UV luminosity 
function derived by Bouwens et al. (2010b, 2010c). The 
UV luminosity function of Bouwens et al. (2010b) was ob- 
tained by using their five T103 -dropout galaxies in HUDF dis- 
cussed in this paper. The best-fit Schechter function parame- 
ters are M* = -19.5 ±0.3, 0* = 1.1 x lO^Mpc" 3 (fixed), and 
a = -1.74 (fixed). Recently, Bouwens et al. (2010c) ob- 
tained a new UV luminosity function for galaxies at z ~ 8 
using 47 F105 or lo98-dropout galaxies found in the addi- 
tional two HUDF09 fields and deep, wide-area Early Release 
Science (ERS). The best-fit Schechter function parameter 
set is M* = -20.28 ± 0.19, 0* = 0.38^ X lO^Mpc" 3 , and 
a = -2.00 ± 0.33. Because the faint-end slope of the UV lu- 
minosity function of Bouwens et al. (2010c) is steeper than 
that of Bouwens et al. (2010b), the contribution of faint galax- 
ies to the ionizing photon production rate density is larger for 
the UV luminosity function of Bouwens et al. (2010c) than 
for Bouwens et al. (2010b). 

The ratio of ionizing photon production rate to UV (rest- 
frame 1500 A) luminosity depends on the metallicity. Since it 
is probable that the metallicity of stars formed in the galax- 
ies at z ~ 8 is lower than the solar value, we investigate 
the following cases; Z = (i.e., the Population III condi- 

5 The selection criteria in Bouwens et al. (2010c) are slightly different 
from those in Bouwens et al. (2010b): they adopted F105 — /125 > 0.45 A 
J125 -#160 < 0.5 for HUDF09, HUDF09-1, and HUDF09-02. Since the ERS 
observations use a different F-band filter (Yq9%), the selection criteria for ERS 
arey 9 8 -/ I2 5 > 1-25 AJ 125 -H m < 0.5. 



Star Formation at Redshift of 8 



5 



tion), Z = 5 x 1(T 4 Z , Z = O.O2Z , and Z = 1Z©. In or- 
der to convert the UV luminosity density to ionizing pho- 
ton production rate density, we adopt the relation between 
them at age of 10 Myi0 for the constant star formation model 
with Salpeter IMF of (M lower , M uppel -) =(1M Q , 100M Q ) calcu- 
lated by Schaerer (2003). In addition, we also adopt another 
Salpeter IMF of (Mi ower , M upper ) =(1M Q , 5OOM ) for Z = 
and Z = 5 x 10~ 4 Z Q . The basic properties used in this analysis 
are given in Appendix B; i.e., the evolution of ionizing pho- 
ton production rate, Qu, rest-frame UV luminosity at 1500 A, 
La(1500), and their ratio, Qu/L\(l500) as a function of age. 
We note here that we do not include any correction for dust 
obscuration. 

The escape fraction required to reionize the universe is sum- 
marized in Table|4]for the various metallicity and upper mass 
cutoff mentioned above. The tabulated values of / esc are 
shown as a function of L* (i.e., stars that contribute to the 
ionizing photons). It is here noted that the stellar mass of an 
L* galaxy is estimated as M star ~3.5x 10 8 M^ for the UV lu- 
minosity function of Bouwens et al. (2010b)Q Note also that 
Labbe et al. (2010) obtained M* - 10 9 M Q for three galaxies 
at z ~ 8 in HUDF from their very deep SpitzerflRAC obser- 
vation. 

It is shown from Table[4]that galaxies at z ~ 8 with the solar 
metallicity can ionize the intergalactic medium if less-massive 
objects with Muv ~ _ 15 also form massive stars. Note, how- 
ever, that such objects have only M star ~ 10 7 M©. Moreover, 
in order to keep the universe ionized, the source-averaged es- 
cape fraction must be as high as / esc ~ 0.9 for the UV lu- 
minosity function of Bouwens et al. (2010b) and ~ 0.7 for 
that of Bouwens et al. (2010c). It seems likely that some 
galaxies have dusty clouds because supernova events could 
enrich heavy elements in the interstellar gas and then cause 
the formation of dust grains. Since dusty clouds make the 
escape fraction smaller (e.g., Yajima et al. 2009), the source- 
averaged escape fraction could be close to unity. Therefore, in 
the case of solar metallicity, the cosmic re-ionization cannot 
be solely achieved by star-forming galaxies at z ~ 8. 

The chemical evolution of theoretical models of galaxies 
suggests that the metallicity of galaxies at z ~ 8 is less than 
0.01Z Q (e.g., Razoumov & Sommer-Larsen 2010). It is then 
more reasonable to investigate the escape fraction in the low- 
metallicity conditions. In the case of Z ~ 0.02Z Q , galaxies 
brighter than 0.1L* could reionize the universe for both UV 
luminosity function. Under the assumption that galaxies are 
the only ionizing sources, their source-averaged escape frac- 
tion would have to be as large as 0.9. On the other hand, 
in the case of Z = 5 x 10~ 4 Z©, we can obtain a more modest 
required escape fraction, / esc ~ 0.5 for ionization by galaxies 
with L > 0.1L* (M stal > 10 8 M Q ). 

Larger escape fractions (up to / esc ~ 0.8) at higher redshift 
(e.g., z ~ 4 - 10) are suggested from theoretical points of view 
(Wise & Cen 2009; Razoumov & Sommer-Larsen 2010). Ya- 

6 It has been often discussed that the ages of the high-z galaxies are esti- 
mated as a few X 100 Myr or shorter (e.g., Finkelstein et al. 2010; Labbe et 
al. 2010). Since the mean stellar mass that ionizing photons is more massive 
than that emit UV continuum, the ratio of Qh/^-uv becomes smaller with in- 
creasing age. If we evaluate the time averaged Qh/^uv between / = and 
100 Myr, the ratio of Qh/^uv becomes 1 /2 and 2/3 of the value at f=10Myr, 
respectively. This trend is found not only for the case of solar metallicity but 
also for those of lower metallicity discussed in this paper. 

7 We use the relation of M st ar/£uv — lO~ 32 M0/(erg s A -1 ) for the con- 
stant star formation model by Schaerer (2003). This relation is valid for age 
with older than 100 Myr. 



jima et al. (2009) investigated the escape fraction for Lya 
emitters (LAEs) and LBGs at z ~ 3 - 7, taking account of both 
the effect of dust extinction and collisional ionization by su- 
perwind shocks. They found / esc ~ 0.07 - 0.47 for LAEs and 
/esc ~ 0.06 - 0. 17 for LBGs. Although the escape fraction de- 
pends on the systemic mass of galaxies, the lower metallicity 
condition (i.e., Z = 5 x 10~ 4 Z Q ) seems to be more reasonable. 

Finally, we comment on the evolution of the escape frac- 
tion from very high redshift (i.e., z ~ 8) to the present day. 
The available observations indicate / esc < 0.01 for z < 1 while 
/esc ~ 0.01 -0.1 forz ~ 1 — 3 (Inoueet al. 2006 and references 
therein, Steidel et al. 2001; Shapley et al. 2006; Siana et al. 
2007; Iwata et al. 2009 and references therein). However, if 
galaxies were the only sources responsible for re-ionization, 
/esc ~ 0.5 would be required. This dramatic evolution can be 
understood if mass assembly processes could be working at 
z ~ 8. In this evolutionary phase, typical masses of galax- 
ies are less massive (i.e., < 10 9 M Q ) than those in the present 
day, the negative feedback from supernovae is more effective 
and thus superwind-driven bubbles help the escape of ioniz- 
ing photons from galaxies at z ~ 8 (Mori & Umemura 2006; 
Wise & Cen 2009). 

6.3. Blue color of UV continuum 
We discuss the origin of very blue color of UV continuum 
found in some galaxies at z ~ 8; i.e., (3 <-3 (fx oc A^). This 
property has been first found in galaxies at z ~ 7 (Bouwens 
et al. 2010a; Finkelstein et al. 2010). Note that local star- 
bursts and star-forming galaxies at z ~ 5 - 6 show f3 ~ -2. 
Then, Bouwens et al. (2010b) also noted that their galaxies at 
z ~ 8 tend to have similar properties (see also Finkelstein et 
al. 2010). Since the Lya emission falls in Jns band for the 
galaxies at z ~ 8, more careful check is necessary to investi- 
gate their UV continuum. In order to investigate this issue, 
we discuss effects of the escape fraction of ionizing contin- 
uum and Lya emission shifted to the /125 window. 

Based on the population synthesis model of Schaerer 
(2003), the UV continuum slope (3 is shallower than -2.5 for 
any metallicity if the nebular continuum emission is taken into 
account. However, in the case of no nebular continuum emis- 
sion, /? becomes steeper than -3 if the metallicity is lower 
than Z < 5 x 1O~ 4 Z0. These properties are shown in Figure|6j 
see also Bouwens et al. (2010a). 

If the escape fraction is small (i.e., / esc ~ 0), the ionizing 
continuum is absorbed by H I gas within the galaxy, causing 
galaxy-scale H II regions. As a result, the nebular continuum 
emission is so strong that the color of UV continuum becomes 
redder; e.g., j3 > —2.5. In the case of Z - 5 x IO^Zq, the 
escape fraction would have to be / esc ~ 0.5 if galaxies with 
L > O.IL* were to keep the universe ionized (see Table |4j>. 
If this is the case, the nebular continuum emission is highly 
suppressed and thus blue colors of UV continuum should be 
observed. 

In Figure|7l we compare the observed colors and the model 
results. For simplicity, the model loci (dotted curve) for 
(3 = -2.5 and -3.5 are shown as a function of redshift from 
z = 7. Here we use the averaged photometry given in Ta- 
ble 3. All galaxies appear to be consistent with /3 = — 2.5 at 
the 2a level, especially if we consider reddening as low as 
Ay ~ 1. However, the observed scatter in this Figure appears 
to be large (see also Bouwens et al. 2010b; Finkelstein et al. 
2010). Since the age of galaxies at z ~ 8 may be from several 
Myr to a few hundred Myr (note that the age of the universe at 
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Z ~ 8 is 600 Myr), the observed color distributions are under- 
stood if the nebular continuum emission is negligibly small 
(see the right panel of Figure[6]i. 

Two galaxies (5 and 6) appear to have much bluer UV con- 
tinuum; i.e., j3 ~ -4. These colors cannot be explained even 
with very young, extremely metal poor stellar populations 
(see again the right panel of Figure|6J. One possibility is that 
strong Lya emission contributes to the /125 band flux. In or- 
der to explain the observed blue J125 -#i60 colors, the equiv- 
alent width of Lya emission must be as large as 250 A. In the 
right panel of Figure [7j we show the loci of model galaxies 
with Lya emission line [EWo(Lya) = 250 A] for j3 = -2.5 and 
(3 = -3.5. It is shown that the observed data of 5 and 6 are 
explained with the effect of Lya emission. Some high-z Lya 
emitters have EWo(Lya) larger than 100 A (e.g., Nagao et al. 
2004, 2005, 2007; Murayama et al. 2007; Shioya et al. 2009). 
We note that a large amount of intergalactic neutral hydrogen 
absorbs a large fraction of the Lya photons, which would re- 
duce the impact of Lya on the colors. If the metallicity of the 
galaxies at z ~ 8 is much smaller than that at z ~ 5 - 6, the 
Lya equivalent width is expected to be much larger than 100 
A; e.g., EWo(Lya) > 500 A (Schaerer 2003). 

Another interesting property found in Figure UJ is that the 
other two galaxies, 7 and 8, show redder colors, correspond- 
ing to [$ ~— 1. If we explain this property by dust extinction, 
the extinction would have to be as large as Ay ~ 2, corre- 
sponding to the rest-frame UV extinction, Auv ~ 5. If this is 
the case, these galaxies are brighter by 5 mag than the esti- 
mated absolute magnitudes; i.e., they would be as bright as 
M uv ~ -24 mag. Yajima et al. (2009) pointed out that the 
escape fraction of ionizing photons is reduced by dust extinc- 
tion in a factor of from a few to 10. For dusty galaxies, the es- 
cape fraction of the ionizing photons becomes small. If these 
two galaxies are highly extincted and have resultantly low es- 
cape fractions, then to achieve the average escape fraction of 
/esc ~ 0.5, the escape fraction of other galaxies at z ~ 8 must 
be nearly unity. If the cosmic re-ionization is completed by 
star-forming galaxies, it is again suggested that lower metal- 
licity conditions are more favorable at z ~ 8. 



7. CONCLUDING REMARKS 

In this paper, we have found that strong emission-line 
galaxies at z ~ 1 .5 - 2 do not work as interlopers for the iden- 
tification of the 20 galaxies at z ~ 8 in the HUDF by the five 
groups (Bouwens et al. 2010b; Bunker et al. 2010; McLure et 
al. 2010; Yan et al. 2010; Finkelstein et al. 2010). Although 
they are free from low-z interlopers, there may be other un- 
expected problems in the procedures in the source detection. 
Taking these points into account, we have constructed a ro- 
bust sample of eight galaxies at z ~ 8 that are all identified 
by at least two groups. Using this sample, we have investi- 
gated the observed bluer color of UV continuum. Based on 
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the two UV luminosity functions obtained by Bouwens et al. 
(2010b, 2010c), we have also investigated the contribution of 
the galaxies at z ~ 8 to the cosmic re-ionization together with 
the escape fraction of ionizing continuum from the galaxies. 
Our results and conclusions are summarized below. 

1 . Based on the UV luminosity functions of Bouwens et 
al. (2010b, 2010c), we estimate the ionizing photon 
production rates for various stellar populations with 
metallicity from Z = (i.e., the Population III condi- 
tion) to Z = Z Q . Comparing these rates to the critical 
values for the completion of cosmic re-ionization, we 
have found for the solar metallicity stars that the escape 
fraction must be close to unity even if we take account 
of ionizing photons from galaxies with ~ 0.01L*. On 
the other hand, in the case of low metallicity condition, 
Z~5 x 10~ 4 Z Q , the modest escape fraction, / esc ~ 0.5, 
from galaxies with L > 0.1L* is required to achieve the 
cosmic re-ionization. We therefore suggest that metal- 
poor conditions, e.g., Z ~ 5 x IO^Zq, are more favor- 
able for star formation in galaxies at z ~ 8. 

2. Such a large escape fraction can be responsible for the 
observed very blue UV continuum color, (3 < -3, for 
some galaxies at z ~ 8. It is also expected that a large 
number of Lya photons can escape from galaxies and 
then contribute to the flux in /125 band. These two fac- 
tors are responsible for the observed blue colors of UV 
continuum. 

In summary, the stellar populations of the galaxies at z ~ 
8 are dominated by significantly low-metallicity stars with 
Z~5x 1O~ 4 Z0. Subsequent explosions of supernovae lead 
to multiple superwinds in the universe at z ~ 8, and then re- 
ionize the universe. All the observational properties of the 
galaxies at z ~ 8 appear to be consistent with this scenario. 

In this paper, we assume that low-metallicity stars resided 
in the galaxies at z ~ 8 would keep the universe ionized, 
there may be some other ionization sources at such high red- 
shift; e.g., pure Population III stars (Choudhury & Ferrara 
2006), kinematic heating (collisional ionization) by supernao- 
vae (Tegmark et al. 1993; see also Miniati et al. 2004), mini 
quasars (Madau et al. 2004), and so on. However, we are 
not able to determine which sources are more dominant for 
the cosmic re-ionizations at z ~ 8. Future observational and 
theoretical investigations will be necessary to resolve this im- 
portant issue. 
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APPENDIX 

OPTICAL DROPOUT CONSTRAINTS TO REJECT THE CONTAMINATION BY STRONG EMISSION-LINE GALAXIES AT Z ~ 1.5 TO 2 

We demonstrate here that the constraint from the lower limit of hi5~Ji25 color is strongest among the upper limit of B435 -/125, 
V(,o(,-J\25, 2775-/125, and zgso-Jns (i.e., non detection of all optical ACS bands). Figure[8]shows the relation between the upper 
limits of ACS bands and 

In Figure|8] we show the three SEDs of a strong emission-line galaxy of J = 29.2 at z = 1.5 for Ay = 0, 1, and 2. The detection 
limits of the ACS imaging of HUDF are also shown in this figure (Bouwens et al. 2010b). As shown here, if such strong 
emission-line galaxies are dust free, they must be detected in all ACS bands. However, in the case of Ay = 2, the predicted 
magnitudes of the galaxy are fainter than the upper limits in B435, V606, and z&so bands while the predicted magnitude of (775 band 
is brighter than the upper limit. Therefore, the lower limit of (775 -/125 color is strongest among the upper limit of £435-/125, 
V606--A25, '775 —^125, and Z850 _ /i25- The magnitude of the model SED in Figure|8](7 = 29.2) may seem fainter than the observed 
galaxy. We note that the limiting magnitude shown in this figure is measured with 0.35" diameter aperture. On the other hand, 
the magnitudes of galaxies at z ~ 8 are the total ones (see the footnote in Section 6.1). 

METALLICITY DEPENDENCE OF THE IONIZING CONTINUUM 

Many authors (e.g., Bunker et al. 2010) have compared the SFR density based on the UV luminosity density with the critical 
SFR density formulated by Madau et al. (1999). In their formulation, the Salpeter IMF with (Mi ower , M uppa )=(lM Q , 100M Q ) 
and the solar metallicity are used. However, it is known that the ionizing photon production rate and the UV luminosity depend 
strongly on the metallicity (e.g., Schaerer 2003). In order to discuss star formation properties in high-z galaxies, it seems better 
to investigate star formation in metal poor environs. 

The upper panel of Figure [9] shows the evolution of ionizing photon production rate (Qh) for a constant star formation rate 
(SFR = lM Q yr _1 ) calculated by Schaerer (2003). It is shown that the ionizing photon production rate becomes larger with 
decreasing metallicity and with increasing upper mass limit. The middle panel of Figure [9] shows the evolution of UV (rest- 
frame 1500 A) luminosity for the same models as the upper panel of Figure [9] For logf < 6.5 the UV luminosity becomes 
fainter with decreasing metallicity. This occurs because the effective temperature becomes hotter with decreasing metallicity. 
Since the ionizing photon production rate density is converted from the UV luminosity density, the dependence of Qh/L\(1500) 
on Z is important to discuss the ionizing photon production rate density in high-z universe. Figure [9] shows the evolution of 
Qn/L\(1500) normalized by that for Z = Z as a function of metallicity. For galaxies with Z > 0.02Z Q , the ratio is smaller than 
1.82. On the other hand, for the case of Population III, the ratio is larger than 2.35 and can reach 4.58 (5.19) for Salpeter IMF 
withM U p = 100M Q (500M Q ). 
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TABLE 1 

Candidate Galaxies at z ~ 8 in the HUDF 



No. 


Bouwens ID 


Bunker ID 


McLure ID 


YanID 


Finkelstein ID 


1 


y-37636015 


YD7 


2079 


z8-Bl 14 


200 


2 


y-37796000 


YD2 


1939 


z8-Bl 17 


213 


3 


y-38135539 


YD3 


1721 


z8-Bl 15 


125 


4 


y-42886345 


YD1 


1765 


zo-KUyz 


819 


5 


y-430862/6 




2841 


r,0 DAOO 

Z5-B088 


653 


6 




YD4 


2487 






7 




YD5 




z8-SD24 




8 






1422 




2055 


n 

y 




1 1JD 








1 n 

1 u 








iO Du4 1 




11 








Z8-B094 




12 








Z8-B087 




13 








Z8-SB27 




14 








Z8-SB30 




15 








Z8-SD05 




16 








Z8-SD02 




17 








Z8-SD52 




18 








Z8-SD15 




19 










800 


20 










640 


21* 


z-44716442" 


zD7 


1107 


Z7-A044 


3022 



Although the galaxy No.21 is selected as z ~ 8 galaxy from the photometric redshift (McLure et al. 2010; Finkelstein et al. 2010), it is selected as z ~ 7 galaxy 
from the dropout method (Oesch et al. 2010; Bunker et al 2010; Yan et al. 2010). In this paper, we categorize this galaxy as z ~ 7 galaxy. 
** Oesch et al. (2010). 



TABLE 2 

Summary of Selection Criteria of z ~ 8 Galaxies. 



Reference Yiqs -Jns vs. /125 -#160 ACS Bands 7i25 (mag) 



Bouwens et al. (2010b) F105 -7i25 > 0.8, J125 -H m < 0.5, No detection (< 2cr) in all bands < 28.6 

and (J l25 -H l60 ) < 0.2 + 0.12(F 105 -7i 25 ) 
Bunker et al. (2010) Y l05 -J l2 s > 1 No detection « 2cr ) in all bands < 28.5 

McLure et al. (2010) photometric redshift 

Yan et al. (2010) Y [05 -J i2s > 0.8 No detection « 2a) in all bands < 29.0 

Finkelstein et al. (2010) photometric redshift 
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TABLE 3 

Summary of Photometric Properties of Our Robust Sample of Galaxies at z ~ 8 



# 


ref. 


Bouwens 


Bunker 


McLure 


Yan 


Finkelstein 


ave. 


1 


/ 

Y-J 
J-H 


> 1.8 
0.0 ±0.3 


28.44 ±0.16 
1.13 ±0.14 
-0.17 ±0.24 


29.01 ±0.38 
0.61 ±0.79 
-0.32 ±0.63 


28 94 + 23 
>0.56 
-0.55 ± 0.43 


28.02 ±0.12 
> 1.03 
-0.27 ±0.1 8 


28.60 ±0.42 
0.87 ±0.54 
-0.26 ±0.27 


2 


J 

Y-J 
J-H 


27.9 ±0.1 
1.6 ±0.9 
-0.3 ±0.2 


27.88 ±0.10 
2.24 ±0.73 
-0.19 ±0.14 


28.33 ±0.23 
> 1.67 
-0.16±0.35 


28.49 ±0.19 
> 1.02 
0.09 ±0.25 


28.70±0.17 

>0.62 
-0.07 ± 0.25 


28.26 ±0.37 
1.92 ±0.74 
-0.13 ±0.18 


3 


j 

Y-J 
J-H 


2« 4 + n 2 
> 2.1 
0.1 ±0.2 


28 (17 + 011 
1.70 ±0.54 
0.00±0.15 


28 41 + 24 
> 1.59 
0.25±0.31 


28 39 + 1 7 
> 1.13 
0.03 ±0.22 


28 fil +0 14 

>0.86 

0.30±0.18 


28 38 + 21 
1.70 ±0.54 
0.14±0.17 


4 


j 

Y-J 
J-H 


28.1 ±0.2 
1.1 ±0.4 
-0.1 ±0.2 


27 70 + 08 
1.19 ±0.25 
-0.28 ±0.14 


28.23 ±0.21 
1.21 ±0.64 
0.12 ±0.29 


2« 45 ±() 18 
0.85 ±0.55 
-0.08 ±0.25 


28 09 ±0 12 
0.73 ±0.30 
-0.02 ±0.17 


28.11 ±0.28 
1.02 ±0.30 
-0.07 ±0.17 


5 


J 

Y-J 
J-H 


28.6 ±0.2 
1.1 ±0.5 
-0.6 ±0.3 




28.98 ±0.37 
> 1.02 
-0.43 ±0.65 


28.86 ±0.22 
>0.67 
-0.53 ±0.40 


28.31 ±0.15 
>0.73 
-0.61 ±0.30 


28.69 ±0.32 
1.10±0.50 
-0.54 ±0.23 


6 


J 

Y-J 
J-H 




28.11 ±0.11 
1.63 ±0.67 
-0.79 ±0.26 


28.59 ±0.28 
1.14±0.81 
-0.59 ±0.52 






28.35 ±0.37 
1.38 ±0.63 
-0.69 ±0.32 


7 


J 

Y-J 
J-H 




28.38±0.14 
1.06 ±0.42 
0.41 ±0.17 




29.04 ±0.30 
>0.50 
0.04 ±0.40 




28.71 ±0.50 
1.06 ±0.42 
0.22 ± 0.34 


8 


J 

Y-J 
J-H 






28.07 ±0.19 
0.83 ±0.42 
0.24 ±0.25 




27.58 ±0.11 
0.83 ±0.32 
0.37 ±0.14 


27.82 ±0.36 
0.83 ±0.26 
0.30 ±0.17 



TABLE 4 

Minimum Escape Fraction to Ionize the Universe 



/esc 



z 

(Z©) 


Supper 

(M ) 


(M & ) 


L > L* 


L>0.1L* 


L > 0.01L* 








Bouwens et al. (2010b) 
M uv < -19.5 


M uv < -17.0 


M uv < -14.5 


1 

0.02 


100 
100 


l 
l 




0.83 


0.85 
0.50 


5 x 10- 4 


100 


l 




0.54 


0.33 





100 


l 




0.31 


0.19 


5 x ur 4 


500 


l 




0.46 


0.28 





500 


l 




0.28 


0.17 








Bouwens et al. (2010c) 












M uv < -20.28 


M uv < -17-78 


M uv < -15.28 


1 


100 


1 






0.69 


0.02 


100 


l 




0.90 


0.41 


5 x nr 4 


100 


l 




0.59 


0.26 





100 


l 




0.34 


0.15 


5 x ur 4 


500 


l 




0.50 


0.23 





500 


l 




0.31 


0.14 
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FIG. 1 . — Example of the SED of a strong emission-line galaxy at z = 1.5. This spectrum corresponds to one model spectrum generated in Section 4; the SED 
with age of 100 Myr. Note that this is also similar to that of I Zw 18 (Izotov et al. 1999). The filter responses of F105W, F125W, and F160W are also shown. 




t (Myr) 

FIG. 2.— Evolution of EWo(H/3), EW ([O II]A3727), and EW ([O III]A5007), for the constant star formation model with Salpeter IMF and Z = O.2Z = 0.004. 
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FIG. 4. — Diagrams of F105 — /125 vs. /125 — H16O (left), and Y105 —Jva vs. 1775 —Jizs (right)- Blue, green, and red curves show the loci of strong emission-line 
galaxies between z = 1.47 and z— 1.81 with Ay = 0, 1, and 2 mag. The observational data of Bouwens et al. (2010b) (a), Bunker et al. (2010) (b), McLure et al. 
(2010) (c), Yan et al. (2010) (d), and Finkelstein et al. (2010) (e) are shown as filled gray circles with error bars. 
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Fig. 4b. — continued. 
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Fig. 4c. — continued. 
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Fig. 4e. — continued. 
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FIG. 5. — Cumulative ionizing photon production rate density (A/i n ) i s shown as a function of Muv for the UV luminosity function of Bouwens et al. (2010b) 
(left panel) and Bouwens et al. (2010c) (right panel). Solid purple, blue, green, and red curves correspond to the models with Z = 0, 5 X IO^Zq, 0.02Zq, 



and IZq, respectively, with the Salpeter IMF of (Mi ower , M U pp a )=(lMQ, 
5 X 10" 4 Z Q , with the Salpeter IMF of (M lower , M u ppe r )=(lM , 500M Q ). 
necessary to keep the universe completely ionized at z = 8: the line for / esl 



IOOMq). Dotted purple and blue curves correspond to the models with Z = and 
Three vertical lines show the critical ionizing photon production rate density that is 
= 1, the dotted line for / esc = 0.5, and the dashed line for / csc = 0.1. 
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FIG. 6. — Evolution of UV continuum slope, 0, 1550 A in the rest frame is shown as a function of age of the system for different metalhcities. These model 
results are taken from the constant star formation models by Schaerer (2003): (left) stellar continuum + nebular emission, and (right) stellar continuum only. The 
meaning of each curve is the same as that in Figure|5] 




FIG. 7. — Comparison between the observed colors and model results in Fios— Jns vs. /125— ffi60 diagram. Left: observed colors are shown together with 
the locus for model galaxies with f} = —2.5 at z > 1 (dotted curve). Right: observed colors are shown together with the loci for model galaxies with /3 = -2.5 
and f} = -3.5 at z. > 1 (dotted curve). The flux contribution of Lya emission with EWo(Lya) = 250A is also taken into account for both models. The reddening 
vector for Ay = 1 mag is at shown the upper-right corner. 




FIG. 8. — Limiting magnitudes for the ACS filter bands in the HUDF (Bouwens et al. 2010b) are shown together with the three SEDs of a star-forming galaxy 
at z = 1.5; Ay = (blue), Ay = 1 (green), and Ay = 2 (red). 
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FIG. 9. — Upper panel: evolution of ionizing photon production rate, Q\\. Middle panel: evolution of UV (rest-frame 1500 A) luminosity, L^(1500). Bottom 
panel: evolution of Qh/£a(1500) normalized by that for Z = Zq. The meaning of each curve is the same as that in Figure|5] 



